clear            // clear all variables
xdel(winsid())   // close all the graphic windows

//------ PARAMETRES ------//
//Hz
Fe = 40000
ordrePow    = 4

//Hz
FcLow   = 3400
FcHigh  = 200

//------ CALCULS ------//
ordre       = (2^ordrePow) -1
disp(ordre)

Fcl = FcLow/Fe
Fch = FcHigh/Fe

[coeffVal,Ampl,freqVal]=wfir("bp",ordre,[Fch Fcl],"ch",[-1 0.03])

coefficients = coeffVal'
Q15Coeff = floor(coefficients*2^15)

plot2d(freqVal*Fe/1000,20*log10(Ampl))

